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ABSTRACT 

Tests of general relativity (GR) are still in their infancy on cosmological scales, but forthcoming experiments 
promise to greatly improve their precision over a wide range of distance scales and redshifts. One such ex¬ 
periment, the Square Kilometre Array (SKA), will caiTy out several wide and deep surveys of resolved and 
unresolved neutral hydrogen (HI) 21cm line-emitting galaxies, mapping a significant fraction of the sky from 
0 < -z ^ 6. I present forecasts for the ability of a suite of possible SKA HI surveys to detect deviations 
from GR by reconstructing the cosmic expansion and growth history. SKA Phase 1 intensity mapping surveys 
can achieve sub-1% measurements of fas out to z « 1, with an SKAl-MID Band 2 survey out to z < 0.6 
able to surpass contemporary spectroscopic galaxy surveys such as DESI and Euclid in terms of constraints on 
modified gravity parameters if challenges such as foreground contamination can be tackled effectively. A more 
futuristic Phase 2 HI survey of ~ 10® spectroscopic galaxy redshifts would be capable of detecting a ^ 2% 
modification of the Poisson equation out to z « 2. 

Subject headings: cosmology; observations — large-scale structure of universe — gravitation 


1. INTRODUCTION 

General Relativity (GR) has been tested to high precision 
over a comparatively narrow range of scales and redshifts. To 
date, the strongest constraints come from Solar System and bi¬ 
nary pulsar experiments, which represents only a small region 
in the space of gravitational potentials and spacetime curva¬ 
ture spanned by gravitational phenomena (Psaltis 2008; Baker 
et al. 2015). The detection of deviations from GR would have 
important ramifications for fundamental physics, so extending 
tests of gravity to other regimes is of great importance. 

This is particularly the case in cosmology, where the dis¬ 
covery of the accelerating expansion of the Universe has high¬ 
lighted the gaps in our understanding of gravitational physics 
on the largest scales. As a result, alternative theories of grav¬ 
ity have proliferated in recent years (Clifton et al. 2012; Joyce 
et al. 2015), with many of them seeking to replace the need 
for a cosmological constant or dark energy in order to ex¬ 
plain the acceleration. Solar System/pulsar experiments have 
been valuable in thinning out the extensive list of contenders, 
as only the theories that are able to almost exactly reduce to 
GR on these scales can survive the observational constraints. 
Most viable modified gravity (MG) theories therefore feature 
‘screening mechanisms’, which cause deviations from GR to 
switch off on small scales (Khoury 2010; Brax 2012), leav¬ 
ing them with significantly different predictions from GR only 
over cosmological distances.* 

Recent searches for observational signatures of MG have 
primarily focused on understanding how alternative gravita¬ 
tional theories affect several key phenomena: cosmic expan¬ 
sion, the growth of large-scale structure, and light propaga¬ 
tion. Distance measurements from the CMB, baryon acous¬ 
tic oscillations, and Type la supernova surveys have been 
used to reconstruct the expansion history with good precision 
(Planck Collaboration 2015b), but GRh-ACDM expansion his- 

p.j .bull @ astro.uio.no 

* Depending on the screening mechanism, deviations from GR can also 
appear in certain non-cosmological regimes, for example inside evolved stars 
(e.g. Chang & Hui 2011) and dwarf galaxies (e.g. Jain & VanderPlas 2011). 


tories are very often contained within the parameter spaces 
of MG models, so many of them cannot be definitively dis¬ 
tinguished from GR based on the background evolution alone 
(e.g. Capozziello et al. 2005; Nojiri & Odintsov 2007; Kunz & 
Sapone 2007). More decisive is the linear growth rate, which 
is increasingly well constrained by redshift-space distortions 
(RSDs) and other peculiar velocity measurements from large- 
scale structure surveys (Samushia et al. 2013; Raccanelli et al. 
2013; Johnson et al. 2015). Precision tests with large weak 
lensing surveys have also been performed (Thomas et al. 
2009; Simpson et al. 2013). 

MG theories also tend to change the effective strength of 
gravity on non-linear scales, fc > 0.1 Mpc“*^, and the modi¬ 
fications can be significant even if linear scales are left rela¬ 
tively unmodified from the GR case (Stabenau & Jain 2006; 
Laszlo & Bean 2008; Brax et al. 2013; Llinares & Mota 

2014) . It is difficult to separate out the contribution of MG 
from other, less exotic processes, however, due to the diffi¬ 
culty of correctly modelling baryonic effects on small scales, 
or model confusion with massive neutrinos (Baldi et al. 2014; 
Barreira et al. 2014). 

Cosmological tests of GR are still in their infancy, and sig¬ 
nificant increases in precision - as promised by forthcom¬ 
ing experiments like Euclid, DESI, and LSST - will be re¬ 
quired before many theories can be ruled in or out with con¬ 
fidence. There are few clues as to where one might expect 
deviations from GR to manifest themselves, and so broaden¬ 
ing the reach of the tests in scale and redshift is also necessary 
- the ultimate goal being to leave MG theories with ‘nowhere 
to hide’. To date, the most stringent cosmological constraints 
have come from linear scales (0.01 ^ ^ 0-1 Mpc“*^), late 

times (the acceleration era, z < 1), and small/intermediate ar¬ 
eas of sky (~ few x 10^ deg®), as well as the CMB (Zuntz 
et al. 2012; Avilez & Skordis 2014; Planck Collaboration 
2015b), but there are many other places one can look. Eor ex¬ 
ample, novel gravitational phenomena can arise at late times 
on scales of order the horizon size, and over wide angular sep¬ 
arations (Hall et al. 2013; Lombriser et al. 2013; Baker & Bull 

2015) . Modifications of expansion and growth could also ap- 
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pear at higher redshifts, due to the non-trivial time evolution 
of (e.g.) an extra scalar field. Adding observations of these 
as-yet unexplored regimes can only improve our chances of 
seeing some kind of anomalous gravitational behaviour. 

The aim of this paper is to investigate how cosmological 
tests of GR can be improved and extended to other regimes 
using a new class of large scale structure survey at radio fre¬ 
quencies. These will use large radio telescope arrays with 
low-noise wideband receivers (Myers et al. 2009) to map 
out the redshift-space matter distribution over a gigantic vol¬ 
ume, providing precision data over significantly wider sur¬ 
vey areas and redshift ranges than has been possible before. 
Our focus will be on the Square Kilometre Array (SKA), a 
planned general-purpose array split over two main sites in 
South Africa and Australia. The first phase of construction, 
due to finish around 2023, will consist of two sub-arrays: 
SKA 1-LOW, a low-frequency aperture array operating at < 
350 MHz; and SKAl-MID, a conventional mid-frequency ar¬ 
ray of 130 dishes equipped with low noise receivers covering 
^ 350 MHz - 14 GHz. A second phase, scheduled for com¬ 
pletion around 2030, will improve the overall sensitivity by a 
factor of ^ 10. 

The SKA will survey large scale structure primarily by de¬ 
tecting the redshifted neutral hydrogen (HI) 21cm emission 
line from a large number of galaxies out to high redshift. 
This can be achieved in two ways: by measuring the 21cm 
line for many individually-detected galaxies (a galaxy red¬ 
shift survey); or by measuring the large-scale fluctuations of 
the integrated 21cm intensity from many unresolved galaxies 
(intensity mapping; IM). The SKA surveys will cover a com¬ 
bined survey volume and redshift range that is significantly 
larger than that of even Euclid and LSST (albeit with vary¬ 
ing, and sometimes lower, sensitivity), and should even begin 
to probe scales of order the horizon size, k ^ T-L. Here, we 
will focus on background observables and linear scales in the 
post-reionisation Universe, 0 < z < 6; the detectability of 
horizon-scale modified gravity effects with the SKA will be 
examined in a forthcoming paper. 

This paper is organised as follows. In Sect. 2, we describe 
a set of large-scale structure observables that can be used to 
test GR, and outline some parametrisations that can be used 
to connect them with MG theories. In Sect. 3 we give an 
overview of the planned SKA HI surveys, and describe our 
Fisher forecasting methodology. We present forecasts for the 
various observables in Sect. 4, and conclude in Sect. 5. 

2. MODIFIED GRAVITY OBSERVABLES 

In this section, we discuss some of the generic effects that 
modifications to GR have on the cosmological background 
and growth of matter fluctuations in several different regimes. 
Our primary focus is on linear sub-horizon scales where, for 
simplicity, we will employ the quasi-static approximation^ 
and ignore wide-angle corrections (e.g. Raccanelli et al. 2013; 
Hall et al. 2013). 

2.1. Background expansion history 

Most modified gravity theories exhibit background dynam¬ 
ics that can deviate from the standard ACDM evolution. A 
simple example is that of models that introduce a new scalar 
degree of freedom; this is allowed to evolve dynamically, and 

^ Time derivatives of metric perturbations and new degrees of freedom 
are assumed to be of order "H, and subdominant to spatial derivatives (e.g. 
Silvestri et al. 2013). 
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Fig. 1.— Redshift coverage of possible SKA HI intensity mapping and 
galaxy surveys (colours), compared with other future galaxy surveys (grey). 
The SPHEREx line shows the redshift range for two galaxy samples. 

can therefore have a non-trivial redshift-dependent equation 
of state (Ratra & Peebles 1988; Caldwell et al. 1998). Much 
depends on the particular structure of the modified theory - 
a variety of choices regarding the shape of scalar held poten¬ 
tials, couplings to the matter sector etc. can be made relatively 
freely, leading to complicated parameter spaces that are often 
difficult to characterise even for a single theory. 

One simple way of modelling the wide variety of possi¬ 
ble modifications to the expansion history is by introducing 
a effective dark energy equation of state parameter, which is 
commonly parametrised as w(a) » uig -f u>a(l — a) (CPL; 
Chevallier & Polarski 2001; Linder 2003). The primary virtue 
of this parametrisation is its simplicity; this functional form 
for ia(a) is rarely an acceptable fit at all redshifts (e.g. Marsh 
et al. 2014), and a variety of potentially important effects (e.g. 
the sound speed for a scalar held) are neglected. Nevertheless, 
forecasts for wq and Wa are a useful way of comparing the rel¬ 
ative performance of different survey methods in reconstruct¬ 
ing the expansion history, and so we adopt this parametrisa¬ 
tion in what follows. Dark energy and modified gravity con¬ 
straints with more general parametrisations have been consid¬ 
ered elsewhere (e.g. Zhao et al. 2015). 

The equation of state, w(a), can be reconstructed from the 
angular diameter distance-redshift relation, Da(z), and the 
expansion rate, 

H (a) = Hq^J Q,m,oO'~^ + ^^DE(a) + (1) 

where the fractional energy density of dark energy in the CPL 
parametrisation is 

5^DE(a) = Hde,o exp [3wa{a - 1)] (2) 

Typically, a number of probes are combined to pin down these 
quantities over a range of redshifts; we will concentrate on ob¬ 
servations of ‘statistical standard rulers’ in the redshift-space 
matter distribution, combined with a high-redshift constraint 
from the primary CMB power spectrum. The most robust 
such standard ruler is the baryon acoustic oscillation (BAG) 
feature in the matter correlation function. If its true physi¬ 
cal scale is known, the observed redshift-space locations of 
the BAG feature in the radial and transverse directions can be 
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used to constrain two independent combinations of Da and H 
as a function of redshift, with results that are remarkably ro¬ 
bust to systematic effects. Other standardisable rulers can be 
constructed from the matter distribution - the overall (‘broad¬ 
band’) shape of the power spectrum contains more informa¬ 
tion than the BAO alone, for instance - but signihcantly more 
careful modelling is required to avoid systematic biases. Op¬ 
timistically, we will focus on the broadband case here. 


should be further restricted by applying viability conditions 
(e.g. Piazza et al. 2014; Perenon et al. 2015). This is not sim¬ 
ple in practise, so to make progress we follow a similar ap¬ 
proach to Simpson et al. (2013); Baker et al. (2014a); Planck 
Collaboration (2015b) and apply the parametrisation 


p(fc,a) 


1 + 


Ode (a) 
Ode.o 



(8) 


2.2. Linear growth and RSDs 

Gravitational infall is the dominant process in structure for¬ 
mation and so, as one might expect, modifying the theory of 
gravity can have wide-ranging effects on the cosmic matter 
distribution. At linear order in perturbations, a relatively gen¬ 
eral approach to incorporating MG effects is to modify the 
growth equation (valid on sub-horizon scales). 

Am+'H^m = fJ,{k,a)AM, (3) 

where Am is the gauge-invariant matter perturbation and 
overdots denote differentiation with respect to conformal time 
(see, e.g.. Baker et al. (2014b) for a derivation). The func¬ 
tion p(fc, a) has been introduced to parametrise time- and 
scale-dependent deviations from GR growth. In the GR limit, 
/r —>■ 1. We have also dehned flM(a) = /E‘^{a), 

where E{a) = El (a)/EI q is the dimensionless expansion rate, 
and EL = aEl. 

The favoured method of constraining the growth history is 
to measure the linear growth rate. 


/(fc,a) 


d log Am 
d log a 


(4) 


as a function of redshift and scale. Inserting (4) into (3), one 
obtains (Baker et al. 2014b) 


dfjk, a) 
d\oga 


-f{k,a) 


/(fc, a) + ^ - - ^Mia)] 


+ A^M{a)fi{k,a). 


(5) 


Equation (5) is readily solved using standard numerical tech¬ 
niques once p(fc, a) and w{a) are given. While one might be 
tempted to try to preserve generality by keeping p(fc, a) com¬ 
pletely arbitrary, large regions of this space of functions are 
likely to map onto contrived or physically implausible mod- 
ihed gravity theories. Instead, it is prudent to impose some 
structure, such as that suggested by the rather general Hom- 
deski class of second-order scalar single-held modihcations 
to GR (Horndeski 1974; Deffayet et al. 2011). These can be 
shown to restrict the form of 7 and ^ to (de Eelice et al. 2011; 
Amendola et al. 2013b; Silvestri et al. 2013) 


l+P3(a)k^ 

( 6 ) 

n ^ l+P3{a)k‘^ 

E(k,a)- 

pAa) +P5[a)k^ 

(7) 


where {pn} are polynomials in scale factor only. This offers 
a strong restriction on the possible form of any scale depen¬ 
dence of the growth rate, but still allows hve essentially ar¬ 
bitrary functions of time. Many choices of the polynomials 
result in theories that are unphysical (e.g. due to instabilities, 
tachyonic modes and so on), so the space of allowed models 


Only two free parameters have been introduced here: sets 

the amplitude of the growth modihcations (A^ 0 for GR), 

and sets the scale at which the growth begins to devi¬ 
ate from scale-independence. The dependence on Ode (a) 
recognises the motivation of modihcations to GR as a possi¬ 
ble explanation of cosmic acceleration, but alternative redshift 
dependences of the modihcations can equally be considered. 
While there has clearly been some loss of generality, Eq. ( 8 ) 
at least attempts to account for redshift- and scale-dependent 
modihcations to growth in a theoretically-motivated way (par¬ 
ticularly the possibility that MG would introduce a new phys¬ 
ical scale, denoted by Baker et al. 2014a), and so we will 
use it as the preferred parametrisation here. Note that other 
methods of parametrising the growth equations have also been 
considered (e.g. Hu & Sawicki 2007; Daniel et al. 2010). 

A number of phenomenological parametrisations are also in 
use throughout the literature. The simplest is (Peebles 1980; 
Linder & Cahn 2007) 

/(a) = HXf(a), ( 9 ) 

where the 7 parameter is known as the growth index, equal 
to 0.55 in ACDMh-GR. Various extensions to this parametri¬ 
sation have been proposed to capture modihcations to the 
growth history in non-standard scenarios, including allowing 
the growth index to depend on redshift or the dark energy 
equation of state (e.g. Steigerwald et al. 2014), or by also 
modifying the overall amplitude of the growth rate (Di Porto 
& Amendola 2008), 

f{a) = nlj{a)[l + T](a)]. (10) 

Eor ease of comparison with previous results, we include 
forecasts for a few variations on these parametrisations as 
well as for Eq. ( 8 ). In particular, we will consider sim¬ 
ple redshift-dependent modihcations to the growth index and 
amplitude of the growth rate, 7 (a) « 7 o + 7 i(l ~ a) 

77(a) « 770 + ??i(l - a). 

The growth rate can be measured in a number of ways, 
mostly by using probes of the peculiar velocity held (since 
V (X f on linear, subhorizon scales) (e.g. Jain & Zhang 
2008; Kosowsky & Bhattacharya 2009; Hellwing et al. 2014; 
Mueller et al. 2015; Johnson et al. 2015). We will concen¬ 
trate on just one method - observations of redshift-space dis¬ 
tortions (RSDs), the anisotropy induced in the galaxy corre¬ 
lation function in redshift-space by coherent peculiar veloci¬ 
ties (Kaiser 1987; Percival et al. 2011). While not as robust 
as BAOs in terms of insensitivity to complicating non-linear 
processes and other systematic effects (Simpson & Peacock 
2010; Macaulay et al. 2013), one can at least model RSDs 
on large scales using linear cosmological perturbation theory 
plus small corrections, and they have already been success¬ 
fully used to test GR (Guzzo et al. 2008; Samushia et al. 2013; 
Beutler et al. 2014). 

Eor small angular separations, the anisotropy induced in the 
redshift-space matter power spectrum is well described by the 
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Kaiser approximation (Kaiser 1987), 

P(k, z)=FRSD(k, z)P(fc, z) (11) 

PRSD(k, z) = (b(z, k) + f{z, (12) 

where p = cos0 is the angle of the wavevector to the 
line of sight, b is the bias of the tracer population with re¬ 
spect to the dark matter, and an exponential term has been 
added to account for the smearing-out of redshift informa¬ 
tion on small scales by incoherent, non-linear peculiar veloc¬ 
ities. The growth rate can be separated from the (generally 
poorly-known) bias by comparing moments of the redshift- 
space matter distribution. There is a degeneracy with the 
normalisation of the power spectrum, as, however, such that 
only the combinations bag and fas (or alternatively bas and 
P = f jb) can be measured directly. This degeneracy can 
be broken using one of the parametric growth models from 
above, combined with a constraint on the normalisation of the 
power spectrum from the CMB. 


3. FORECASTING FOR SKA HI SURVEYS 

In this section, we describe the HI galaxy redshift and in¬ 
tensity mapping surveys that will be performed by the SKA, 
and outline the formalism used to forecast constraints on the 
expansion and growth rates and their parametrisations. 

Galaxy redshift surveys are a tried and tested technique in 
the optical and near-infrared, where they have already been 
used to measure the BAO and RSDs to high precision at 
z < 1. Their survey speed scales poorly with increasing vol¬ 
ume, however, as they rely on making time-consuming high- 
SNR detections of individual objects. Intensity mapping re¬ 
laxes the requirement for individual detections, thus promis¬ 
ing dramatically improved survey speeds. It is so far a rel¬ 
atively untested method however, with the first detection of 
the cosmological fluctuations having been made only recently 
with this technique (Masui et al. 2013). Nevertheless, a host 
of IM experiments are planned to demonstrate the feasibil¬ 
ity of the method over the next few years (Bull et al. 2015b), 
paving the way for the surveys with Phase 1 of the SKA that 
we consider here. We will consider the two methods on an 
equal footing. 


3.1. Fisher forecasting formalism 

Fisher forecasting is a simple, computationally inexpen¬ 
sive way of predicting the constraints on a set of parameters 
that should be achieved by a given experimental configura¬ 
tion. While clearly approximate and idealised - it assumes 
Gaussianity and neglects systematic biases - Fisher forecast¬ 
ing is nevertheless a reliable way of understanding the relative 
performance of different experiments and getting a handle on 
correlations between parameters. 

To proceed, one must first define fiducial models for the ex¬ 
pected signal and noise for a set of observations, as a function 
of the parameters of interest. Models for galaxy redshift and 
intensity mapping surveys are constructed in subsequent sec¬ 
tions, based on the formalism developed in Bull et al. (2015b). 
The Fisher matrix for a set of parameters {0} can be written 
as (Seo & Eisenstein 2007) 


F- ■ — 


(Pk 


Kff(k) 


a log (7^9 log c® 
36, d9~ 


(13) 



Fig. 2.— The sensitivity of the ‘Rebaselined’ and ‘Alternative’ designs for 
the SKAl-MID receivers, shown as total Aeff/Tgys curves for the sub-array. 
The MeerKAT bands (orange lines) are shown for comparison. The yellow 
shaded regions denote the SKAl Baseline specification sensitivity (Dewd- 
ney et al. 2013) for the original assumption of 190 dishes (upper limit), and 
corrected for the post-rebaselining figure of 130 dishes (lower limit). 


The effective volume. 


Ueff(k) = /skyU 


C^(k,z) 

CS(k, z) + C'^(k,z) 


(14) 


is a weighting that accounts for the varying sensitivity of an 
experiment to different Fourier modes (e.g. due to instru¬ 
mental beam effects or cosmic variance), and depends on the 
physical volume of the redshift bin, Vi = [^““(fiWdzjiiz, 
and the fraction of the sky covered by the survey, /sky = 
5'area/47r. The C''^(k, z) and C^jk, z) terms are the signal 
and noise covariance respectively, and will be explicitly de¬ 
fined in subsequent sections. After calculating the Fisher ma¬ 
trix, one can invert it to get an estimate of the expected co- 
variance between the parameters, F~^ « Cov({0}). 


3.2. SKA array configurations 

The design of the SKA has not yet been finalised, so there 
is some freedom in what to assume for the instrumental spec¬ 
ifications. The most complete specification for Phase 1 of 
the SKA is currently the ‘baseline’ design of Dewdney et al. 
(2013). This originally called for a three-array system, con¬ 
sisting of two dish arrays with mid-frequency receivers cover¬ 
ing multiple bands (SKAl-MID and SUR) and a single low- 
frequency aperture array (LOW). This specification has now 
been updated following a ‘rebaselining’ procedure however 
(McPherson 2015), which has removed the SUR array, halved 
the number of receiving stations of LOW, and reduced the 
number of MID dishes by 30% in order to meet a cost cap. 

Following this, options for redefining the available MID 
and LOW bands are also being considered, with the inten¬ 
tion of better-aligning the frequency coverage with recently- 
selected ‘key science’ goals (Braun & Cornwell 2014). This 
is also an opportunity to improve the overall sensitivity of the 
array: by reducing the bandwidth of Band 1 of MID, simpler 
low-noise receivers (with a maximum:minimum frequency ra¬ 
tio of 1.85:1 or less) can be used instead of the more complex, 
higher-noise wideband (3:1 ratio) receivers of the current de¬ 
sign, which also suffer from degraded performance at lower 
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Telescope array 


^recv 

Aeff 

-^dish 

Ng X Ng 

FOV 

^ref 

^min 

^max 

■^min 

^max 

'S'area IM 

Sarea GS 



[K] 

[m2] 

[m] 


[deg2] 

[MHz] 

[MHz] 

[MHz] 



[deg2] 

[deg2] 

SKAl-LOW 

Rebaselined 

40 

925 

35 

455 X 3 

28 

110 

50 

350 

3.06 

27.4 

1,000 



Alternative 

40 

925 

35 

455 X 3 

28 

no 

50 

500 

1.84 

27.4 



SKAl-MID B1 

Rebaselined 

23 

133 

15 

130 X 1 

1.78 

700 

350 

1050 

0.35 

3.06 




Alternative 

12 

133 

15 

130 X 1 

1.78 

700 

450 

825 

0.72 

2.16 

25,000 

— 

MeerKAT 

UHF-band 

23 

115 

13.5 

64 X 1 

2.20 

700 

580 

1015 

0.40 

1.45 



SKAl-MID B2 

Rebaselined 

15.5 

150 

15 

130 X 1 

0.87 

1000 

950 

1760 

0.00 

0.49 




Alternative 

15.5 

150 

15 

130 X 1 

0.87 

1000 

795 

1470 

0.00 

0.79 

25,000 

5,000 

MeerKAT 

L-band 

30 

122 

13.5 

64 X 1 

1.08 

1000 

900 

1670 

0.00 

0.58 



SKA2 

15 

6 

3.1 

70000 X 1 

30 

1000 

470 

1290 

0.10 

2.00 

— 

30,000 


TABLE 1 

Representative instrumental parameters eor various assumed SKA Phase 1 and 2 array configurations and bands. Agg is the 

EFEECTIVE AREA PER DISH OR RECEIVING ELEMENT, AND THE FOV IS CALCULATED AT THE REFERENCE FREQUENCY, U^eC- L'min,max ARE THE BAND 
EDGES, AND Zmin.max THE CORRESPONDING MINIMUM/MAXIMUM ACCESSIBLE REDSHIFTS. REPRESENTATIVE SURVEY AREAS FOR INTENSITY 
MAPPING (IM) AND GALAXY SURVEYS (GS) ARE GIVEN IN THE EINAL TWO COLUMNS. OTHER QUANTITIES ARE DEEINED IN THE TEXT. 


frequencies. Fig. 2 shows a measure of sensitivity - the total 
effective area (^eff) of the array divided by the system temper¬ 
ature {Tsys) - for the original baseline specihcation, the most 
up-to-date sensitivity estimates following rebaselining (SKA 
Organisation 2015), and a proposed ‘alternative’ Band 1 and 
2 conhguration with simpler, reduced-bandwidth receivers (J. 
Jonas & M. Santos, priv. comm.). The alternative conhgura¬ 
tion almost doubles Aes/Tsys at 800 MHz, mitigating the cuts 
from rebaselining. The maximum frequency of LOW could 
be increased to 500 MHz to compensate for the resulting loss 
of frequency coverage of Band 1, without affecting the sensi¬ 
tivity of LOW in its design range of 50 < iz < 350 MHz (K. 
Zarb-Adami, priv. comm.). 

To rehect the Huidity of the design, we therefore chose two 
representative conhgurations for the SKA 1-LOW and MID 
arrays: a ‘Rebaselined’ design based on the recommendations 
of McPherson (2015) and using the updated sensitivity esti¬ 
mates of SKA Organisation (2015); and an ‘Alternative’, also 
based on the rebaselined specihcation, but with a different set 
of frequency bands and reduced system temperature to rehect 
the expected improvement in sensitivity if simpler receivers 
can be used. In our forecasts for SKA 1-MID we also include 
the 64 MeerKAT^ dishes, which signihcantly increases the 
sensitivity over part of each band (see Fig. 2). The specihca- 
tions are summarised in Table 1 . 

The design for Phase 2 of the SKA is yet to be formally 
dehned, and only notional specihcations currently exist.^^ The 
expectation is that it will have approximately 10 x the sensi¬ 
tivity of SKAl above 350 MHz, and will be capable of (spec¬ 
troscopically) detecting ~ 10® HI galaxies. In this paper, we 
assume a dense aperture array design (c.f. Torchinsky et al. 
2015) with the frequency range used in Yahya et al. (2015), 
and choose the collecting area and number of receivers such 
that it would detect ~ 10® HI galaxies above a lOcr threshold 
(see Sect. 3.3). The resulting specihcation is listed in Table 1. 
Note that this actually leads to a total A^s/Tsys « 15x that 
of SKAl-MID H- MeerKAT at 1 GHz. 

^ http : //public . ska . ac . za/meerkat 
http : / / astronomers . skatelescope . org/ska2/ 


3.3. Galaxy redshift surveys 

Galaxy surveys detect and then measure the positions and 
redshifts of many individual galaxies, which are assumed to 
Poisson sample the underlying dark matter density held. The 
measured galaxy correlation function can then be used to infer 
the matter power spectrum in redshift-space, which contains 
a host of information about the growth of structure and other 
cosmological quantities, as discussed above. 

The signal covariance for a galaxy redshift survey is simply 
the redshift-space power spectrum, 

C^ = P(k,z)=FRSD(k,z)P(fc), (15) 


where Prsd was dehned in Eq. (12), and we marginalise over 
the non-linear velocity dispersion scale, unl ~ 7 Mpc, as 
a nuisance parameter. The noise covariance is modelled as 
shot noise, = l/n{z). In the limit that n{z) —/ oo, 
14ff fskyVi, and the Fisher errors are limited only by the h- 
nite number of Fourier modes in the survey volume (i.e. sam¬ 
ple variance). We set fcmax = 0.14 Mpc“^(l -f 
corresponding to a non-linear cutoff on small scales (Smith 
et al. 2003). 

The number density and bias functions depend on the type 
of galaxy being targeted by the survey, as well as the sensi¬ 
tivity of the survey instrument. The rms hux sensitivity for 
dual-polarisation receivers is given by 




2fcBTs, 




(16) 


where tp is the integration time per pointing, N^mh and A^ff 
are the number and effective area of the collecting elements, 
Siy is the channel bandwidth, and the total system temper¬ 
ature is Tsys = Trecv + Tgky, where the hrst term is the 
total contribution of the receiver system to the noise, and 
the second is due to background emission, Tsky « 60 K x 
(iz/300MHz)“®-®. For all array conhgurations, we assume a 
per-element collecting area of A^s « e7r(£)dish/2)®, where 
Pdish is the dish/station diameter and e « 0.7 — 0.9 is a 
typical aperture efficiency. The integration time per point- 
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SKAl-MID B2 + MeerKAT 


■^min 

■^max 

n{z) [Mpc 

b{z) 

‘S'rms [/^Jy] 


0.0 

0.1 

2.73 X 10-2 

0.657 

117.9 

< 

0.1 

0.2 

4.93 X 10-3 

0.714 

109.6 


0.2 

0.3 

9.49 X 10-"^ 

0.789 

102.9 

s 

OJ 

0.3 

0.4 

2.23 X 10-^^ 

0.876 

97.5 


0.4 

0.5 

6.44 X 10-5 

0.966 

93.1 


0.5 

0.6 

1.70 X 10-5 

1.059 

89.4 

< 

0.6 

0.7 

1.72 X 10-5 

1.208 

106.2 


TABLE 2 

Binned number density and bias oe HI galaxies, and 
CORRESPONDING ELUX SENSITIVITY, EOR SKAl-MID BAND 2, 

ASSUMING A 5,000 DEG^, 10,000 HOUR SURVEY. THE DETECTION 
THRESHOLD IS 5(7. 

ing can be rewritten in terms of the total survey time, tp = 
ttot FOV/S'area, where the instantaneous field of view is 

FOV« f(1.3A/Z?dish)". 

The number density and bias of HI galaxies as a func¬ 
tion of redshift and flux sensitivity were calculated by Yahya 
et al. (2015) using the S-cubed simulations (Obreschkow et al. 
2009), which are based on the Millennium dark matter only 
simulation (Springel et al. 2005). A frequency resolution of 
10 kHz and detection threshold of lOcr was assumed in these 
calculations, and the bias was calculated using the halo model 
with a mass function measured from the simulations. We cal¬ 
culate the number density and bias for the array configura¬ 
tions in Table 1 by interpolating the results of Yahya et al. 
(2015) for 10,000 hour surveys over either 5,000 deg^ (SKAl) 
or 30,000 deg^ (SKA2). The dependence of survey perfor¬ 
mance on the assumed survey time and area is discussed in 
Appendix A, and fitting functions for the number density and 
bias are given in Appendix B. To boost the number density 
of detected sources, we also revise the detection threshold for 
SKAl to 5(7, which is relatively low (c.f. Abdalla & Rawlings 
2005). Only Band 2 of MID is usable; the number densities 
obtained with Band 1 and LOW would be too small for a suffi¬ 
ciently large survey area and reasonable survey time. We have 
also assumed a source detection efficiency of 100%, but note 
that effects like source confusion will reduce this (Jones et al. 
2015). The resulting values of n{z) and b{z) in bins of width 
Az = 0.1 are listed in Table 2 (SKAl) and Table 3 (SKA2). 

3.4. Intensity mapping surveys 

The process of cataloguing millions of galaxies over a large 
volume is extremely time consuming, as each source must be 
detected at a sufficiently high signal-to-noise ratio (SNR) to 
ensure that it is not just a statistical fluctuation. An alter¬ 
native approach is to make low-resolution maps of the in¬ 
tegrated emission from many galaxies, in a process known 
as intensity mapping (Madau et al. 1997; Bharadwaj et al. 
2001; Battye et al. 2004; Loeb & Wyithe 2008; Chang et al. 
2008). Fluctuations in the observed intensity of redshifted 
HI emission follow fluctuations in the underlying matter den¬ 
sity field as traced by the HI emitting galaxies, allowing the 
density field to be reconstructed on sufficiently large scales 
from intensity maps. Redshift information is included au¬ 
tomatically as the target is an emission line (i.e. the fre¬ 
quency of observation maps directly to the redshift of emis¬ 
sion, z^obs = *^ 2 icm/(l + z)), and one only need integrate 
down to SNR « 1 to (statistically) detect the cosmological 
signal (c.f. measurement of CMB anisotropies). Intensity 



SKA2 


■^max 

n{z) [Mpc-3] 

b(z) 

‘S'rms [/^Jy] 

0.1 

0.2 

6.20 X 10-2 

0.623 

6.1 

0.2 

0.3 

3.63 X 10-2 

0.674 

5.7 

0.3 

0.4 

2.16 X 10-2 

0.730 

5.4 

0.4 

0.5 

1.31 X 10-2 

0.790 

5.2 

0.5 

0.6 

8.07 X 10-3 

0.854 

5.0 

0.6 

0.7 

5.11 X 10-3 

0.922 

4.8 

0.7 

0.8 

3.27 X 10-3 

0.996 

4.7 

0.8 

0.9 

2.11 X 10-3 

1.076 

4.6 

0.9 

1.0 

1.36 X 10-3 

1.163 

4.6 

1.0 

1.1 

8.70 X 10-"^ 

1.257 

4.5 

1.1 

1.2 

5.56 X 10-4 

1.360 

4.5 

1.2 

1.3 

3.53 X 10-4 

1.472 

4.5 

1.3 

1.4 

2.22 X 10-4 

1.594 

4.5 

1.4 

1.5 

1.39 X 10-4 

1.726 

4.5 

1.5 

1.6 

8.55 X 10-5 

1.870 

4.5 

1.6 

1.7 

5.20 X 10-5 

2.027 

4.5 

1.7 

1.8 

3.12 X 10-5 

2.198 

4.6 

1.8 

1.9 

1.83 X 10-5 

2.385 

4.6 

1.9 

2.0 

1.05 X 10-5 

2.588 

4.7 


TABLE 3 

Number density and bias oe HI galaxies eor SKA2, assuming a 
30,000 deg^ survey eor 10,000 hours, and a lOo- threshold. 


mapping surveys are therefore more rapid, making it possi¬ 
ble to probe larger volumes in a reasonable amount of time. 

The signal covariance for the HI brightness temperature 
fluctuations measured by an intensity mapping survey is 

C'''(k, z) = T^FnsuiK z)Pik), (17) 


where the only difference with Eq. (15) is a factor of the HI 
brightness temperature squared, T^{z). We use the bias and 
brightness temperature models from Bull et al. (2015b). The 
noise covariance for an IM survey can be written as a scale- 
dependent effective ‘number density’ (Bull et al. 2015b), 


l^IM(k., z) —Auish ^hea,m kIpo\ 


X 



^ ^tot 
•^area^HI 


BlBu. 


( 18 ) 


This expression is based on the ideal radiometer equation. 
The factors of B± and Ry are window functions in Fourier 
space due to the angular beam and hnite frequency channel 
bandwidth, and Adish, IVbeam, and Apoi are factors that ac¬ 
count for the gain in survey speed for a telescope array with 
multiple dishes/beams/polarisation channels. The noise also 
depends on the survey area, 5'areaj and the total observation 
time for the survey, ftot- See Table 1 for the values of all 
relevant parameters for the various SKA configurations. 

Eq. (18) depends on the type of receiver system through the 
B± term. We have assumed that the SKAl-MID IM surveys 
will be performed in autocorrelation mode; in interferometer 
mode, MID has too small a field of view and an insufficient 
density of short baselines to be sensitive to the large scales we 
are interested in at redshifts of z < 1.4 (Bull et al. 2015b). 
Its interferometric resolution is better matched to (e.g.) the 
BAO scale at higher redshifts, however. Autocorrelation ob- 
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servations are sensitive to all scales between approximately 
the size of the beam and the survey area, so are better matched 
to the BAG scale at lower redshift for MID. They suffer from 
correlated (1//) noise and ground pickup however, which can 
significantly increase the difficulty of recovering the cosmo¬ 
logical signal (although foreground cleaning and appropri¬ 
ate scanning strategies can help to remove these effects; see 
Bigot-Sazy et al. 2015). SKAl-LOW is an aperture array, 
which operates as an interferometer and has a configuration 
that is better-matched to the BAG scale at high z. Noise and 
beam expressions for the various types of receiver, includ¬ 
ing baseline density distributions for the interferometers, are 
given in Santos et al. (2015) and Bull et al. (2015b). 

Another important systematic effect is the presence of fore¬ 
ground contamination. Galactic synchrotron and other fore¬ 
grounds are around 5-6 orders of magnitude brighter than the 
cosmological HI signal, and so must be removed with a high 
level of efficiency. Most foregrounds should be spectrally 
smooth, making it possible to subtract them using polyno¬ 
mial fitting. Principal Component Analysis, or similar (Wang 
et al. 2006; Morales et al. 2006; Liu & Tegmark 2011; Petro- 
vic & Gh 2011; Alonso et al. 2015). The frequency depen¬ 
dence of the beam response can hinder this process, however, 
with interferometers in particular susceptible to the genera¬ 
tion of non-smooth foreground signals due to chromatic/wide- 
held effects (the ‘foreground wedge’: Bowman et al. 2009; 
Liu et al. 2009; Datta et al. 2010; Morales et al. 2012; Thya- 
garajan et al. 2015; Seo & Hirata 2015). Gther effects, such 
as atmospheric noise (Bigot-Sazy et al. 2015), ionospheric 
distortions, and radio-frequency interference (Alonso et al. 
2014), can also be problematic at low and high frequencies 
respectively. Nevertheless, recent simulation work has shown 
that existing foreground removal methods can recover the true 
HI power spectrum to within 5%, although over-subtraction 
of the HI signal biases the recovered spectrum in a scale- 
dependent way (Wolz et al. 2014; Alonso et al. 2014; Bigot- 
Sazy et al. 2015; Glivari et al. 2015). This is problematic if 
the aim is to use the broadband shape of the power spectrum 
for cosmology; recovery of the BAG scale is not biased by 
this effect, however (Wolz et al. 2014; Alonso et al. 2014). 

We assume ftot = 10'* hours for all IM surveys. The survey 
area is taken to be 25,000 deg^ for MID, and 1,000 deg^ for 
LGW. The MID results do not depend strongly on the assumed 
survey area for S'area ^ 5, 000 deg^; see Appendix A. 

3.5. Prior information 

It is useful to include prior information from other sources 
in the forecasts, e.g. in order to break degeneracies. By the 
time of the hrst SKA HI surveys in the early 2020’s, a large 
amount of precision data from various sources will already be 
available. Rather than trying to forecast for the entire state of 
observational cosmology at that time, we take a more conser¬ 
vative approach and restrict the prior information to just two 
sources: the CMB angular power spectrum from Planck, and 
galaxy clustering information from BGSS, which anchor the 
constraints at high- and low-redshift respectively. Lor Planck, 
we use the DETL Lisher matrix prior, calculated assuming full 
polarisation, 80% sky coverage, and 3 frequency bands free of 
foreground contamination (Albrecht et al. 2009). Lor BGSS, 
we take the binned number density and bias values from Lont- 
Ribera et al. (2014) for a 10,000 deg^ survey and perform 
our own forecasts using the procedure outlined in Section 3.3. 
The priors are applied by adding the Planck and BGSS Lisher 
matrices to the Lisher matrix for a given SKA survey. 


3.6. Parameters used in the forecasts 

In all cases we forecast for the parameters 

{DAiz),H{z),f(Ts{z),bas{z),crnL}, 

where Da is the angular diameter distance and H is the ex¬ 
pansion rate. The hrst 4 parameters are assumed to be free in 
each redshift bin, and the non-linear velocity dispersion, (Tnl, 
is marginalised as a nuisance parameter. This set of param¬ 
eters can be viewed as “model-independent”, as we have not 
assumed parametric functional forms for any of the hrst four 
functions. No priors are applied to this set. 

Lor the forecasts where a growth parametrisation is as¬ 
sumed, we project from the base parameters to 

{h, Gde, Wa, ris, erg, b{z), ctnl} + {^'mg}, 

where {0 mg} growth parameters from one of the 
parametrisations discussed in Sect. 2.2. Both the BGSS 
and Planck priors are applied to this set, and the bias is 
marginalised over as a free parameter in each redshift bin. 
With a parametric model chosen for f{z), the degeneracy be¬ 
tween the bias and normalisation of the power spectrum that 
occurs in the RSD term is now broken, as the functional form 
of crg(z) can be calculated from the growth model, and its nor¬ 
malisation is set by the measured CMB normalisation.^ Ex¬ 
plicitly, we write crg(z) = crgGjz), where erg = crg(z = 0) 
is now a separate parameter and D{z) is the linear growth 
factor. Since D{z) depends on / through the dehnition 
f = d log D/d log a, its derivatives with respect to the growth 
parameters must also be taken into account. We do this by 
projecting /erg into erg and the growth (and CPL) parameters. 
Derivatives of &erg with respect to these parameters are ne¬ 
glected in the projection, as this quantity would not be used 
to constrain {0 mg} in n realistic analysis. We continue to 
marginalise over the bias by projecting bag —t b only. 

In all forecasts, we assume that information about the 
full shape of the power spectrum can be reliably recovered. 



■.I 

0 1 2 3 4 5 6 


z 

Fig. 3.— Derivatives of f{z), E = H{z)/H q, and Dp^{z) with respect to 
the modified growth and equation of state pai'ameters. The Da curves have 
been rescaled by a factor of 2iTo/c. 

^ For IM surveys, these quantities are also degenerate with Ti){z). We 
assume that this is already known, and can be fixed in our analysis. 
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2 : 2 : 

Fig. 4.— Forecast constraints on H{z'] and Dp^(z) for various galaxy surveys (GS) and intensity mapping surveys (IM). Dashed lines show the Alternative 
array configurations. The Ho survey, BOSS and HETDEX specifications were taken from Amendola et al. (2013a) and Font-Ribera et al. (2014) respectively. 


Whether this is possible will depend on having sufficient con¬ 
trol over scale-dependent systematic effects (e.g. see Sect. 
3.4). A more conservative approach would be to use only 
RSDs and BAOs, as the latter are very robust to systematics 
(Seo et al. 2010; Mehta et al. 201 1). 

Finally, we use the Planck Collaboration (2015a) 
flat ACDM cosmology as our fiducial model, 
{nM,^lb,Kns,as) = (0.316,0.049,0.67,0.962,0.834), 

unless otherwise specified. 

4. RESULTS 

In this section we present forecast constraints on the back¬ 
ground expansion rate, linear growth rate, and associated ex¬ 
pansion and growth parametrisations. 

To frame the discussion, it is hrst useful to divide the post¬ 
reionisation Universe into four approximate redshift regimes, 
each of which has something different to say about possi¬ 
ble deviations from GRh-ACDM. To help illustrate the dif¬ 
ferences between the regimes. Fig. 3 shows the derivatives of 
f{z), H{z), and Da(z) with respect to the expansion/growth 
parameters as a function of redshift for the hducial cosmol¬ 
ogy. For a given measurement precision, measurements of the 
growth and expansion rates should be more sensitive to these 
parameters where the derivatives are largest, although de¬ 
generacies with other parameters strongly influence the con¬ 
straints as well. The four regimes are as follows; 

Low redshifts (0 < 2 < 0.7) — The expansion becomes dark 
energy-dominated and the linear growth rate deviates appre¬ 
ciably from unity here, making this the most obvious regime 
in which to look for modifications to GR that are motivated as 
an explanation for cosmic acceleration. Fig. 3 suggests that 
this redshift range should be the most sensitive to the growth 
index parameters, 70 and 71 . Existing surveys have already 
made ~few percent-level measurements of Da{z) and H{z) 
(~ 10 % for fas) at these redshifts which, in combination with 
the CMB, constitute some of the most stringent cosmological 
constraints on GR to date (Planck Collaboration 2015b). 

Dark energy transition (0.7 z < 2) — This is where dark 
energy does not yet dominate the expansion but has started to 
become dynamically important, making it a suitable regime 
to look for evolution of the equation of state and growth in¬ 


dex. This is backed up by Fig. 3, which shows that the deriva¬ 
tives peak in this range for many of the parameters. Precision 
measurements of expansion and growth have only reached the 
lower part of this range so far, but it is an important focus for 
a number of forthcoming surveys like Euclid and DESI. 

Intermediate redshifts (2 < 2 < 3.5) — The expansion is 
matter-dominated in this range - dark energy is a subdomi¬ 
nant (< 10%) contribution to the net energy density of the 
Universe, and the growth rate approaches unity in GR. Nev¬ 
ertheless, the high-redshift behaviour of the equation of state 
can be probed here, w{a) —)■ u>o + Wa, as can a rescaling 
of the growth rate, 1 - 1 - 77 ( 0 ) —>■ 1 -f 770 + t?i- There is al¬ 
ready a ~ few percent distance constraint at 2 = 2.4 from the 
BOSS Lyman-a BAO measurement (Slosar et al. 2013), and 
the planned HETDEX survey® should provide similar preci¬ 
sion on expansion and growth over most of this range. 

Post-reionisation (3.5 < 2 ; < 6)— This is physically similar 
to the previous redshift range except that the galaxy popula¬ 
tion is younger. While some MG/dark energy theories may 
modify the expansion rate at these high redshifts (Chiba et al. 
2013; Tsujikawa 2013), this regime is more likely to be use¬ 
ful as a consistency check on the lower redshift constraints, 
e.g. as a way of resolving degeneracies between massive neu¬ 
trinos and modified gravity (Baldi et al. 2014; Barreira et al. 
2014; Villaescusa-Navarro et al. 2015). There are no existing 
expansion/growth constraints at these redshifts. 

The specihcations of the surveys that will cover the vari¬ 
ous regimes can be found in Sects. 3.2-3.4 (summarised in 
Table 1), and their redshift coverage is illustrated in Fig. 1. 

Note that forecast constraints on Da{z), H{z), the dark 
energy equation of state, and fas, have been presented be¬ 
fore for several of the SKA surveys listed here - for ex¬ 
ample, in Bull et al. (2015a), Raccanelli et al. (2015), and 
Yahya et al. (2015). The key difference in this work is that 
the specihcations have been substantially updated following 
the SKA rebaselining procedure (Sect. 3.2). Specihcally, the 
previous works used the baseline specihcations for MID and 
LOW, which have seen reductions in collecting area of ap¬ 
proximately 30% and 50% respectively. 

® http : //hetdex . org/ 
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BOSS 
Ha survey 
MID B1 + MK (Alt.)" 
MID B2 + MK (Alt.) 
LOW (Alt.) 

SKA2 



0.70 
0.65 
0.60 
f- 0.55 
0.50 
0.45 
0.40 


T- 1 -r 



J_ I _I_ L 


-1.4 -1.2 -1.0 -0.8 -0.6 

Wq 


Fig. 5.— Joint constraints on wq, Wa, and 7 for the 7 ( 2 ) = 70 growth parametrisation. All of the ellipses include the BOSS and Planck Fisher matrix priors. 
All other parameters have been marginalised. 


4.1. Background expansion 

Fig. 4 shows the forecast constraints on the expansion rate, 
H{z), and angular diameter distance, Da{z). At low red- 
shifts, an HI galaxy survey with SKAl-MID Band 2 will per¬ 
form worse than current optical surveys owing to its low sur¬ 
vey volume, as was discussed in Yahya et al. (2015). A 25,000 
deg^ intensity mapping survey in autocorrelation mode with 
SKAl-MlD Band 2 will provide highly competitive ^ 1% 
constraints out to z « 0.7 however, assuming that fore¬ 
grounds and other systematics can be handled without causing 
too large a loss in effective sensitivity. 

Similar precision on H{z) could be achieved out to z « 2 
with an equivalent autocorrelation-mode IM survey on MID 
Band 1, although this is less competitive because Ha surveys 
like Euclid will provide signihcantly better (~ 0.4%) con¬ 
straints in the same redshift range in the same timeframe. The 
MID Band 1 constraints on Da degrade more rapidly with 
redshift due to the resolution effects discussed in Sect. 3.4, 
but remain at the few-percent level. Note that the Alternative 
SKAl-MID conhguration would result in a small improve¬ 
ment in H and Da constraints over the Rebaselined specifi¬ 
cation, although the loss of redshift bins at z < 0.7 would 
decrease the MID Band 1 dark energy figure of merit. 

More promising is an SKA2 HI galaxy survey, which will 
measure H{z) with sub-percent precision over 30,000 deg^ 
from z « 0.1 — 2.0, reaching the 0.3% level around z « 1 
where sensitivity to the dark energy equation of state param¬ 
eters peaks (according to Fig. 3). This supports the classih- 
cation of SKA2 as a key Stage IV dark energy survey by the 
Dark Energy Task Force (DETF; Albrecht et al. 2006) if the 
hnal specihcations are comparable to what we assumed here. 

At higher redshifts (z > 2), IM surveys on SKAl-MID 
Band 1 yield 1 — 3% precision on H{z) for the Rebaselined 
conhguration, at least equalling the performance of HETDEX 
in the same range (albeit with a lower maximum redshift). 
The Da{z) constraint is worse than for HETDEX however, 
and the Alternative MID Band 1 conhguration would not be 
able to reach beyond z = 2 at all. SKA 1-LOW is less sensi¬ 
tive, producing only 5 — 6 % constraints on H{z), but covers 
a wide (and as-yet unprobed) redshift range, so could be use¬ 
ful in constraining MG/DE models that deviate strongly from 


GRh-ACDM only at high redshift. Recall that the performance 
of IM surveys is contingent on the efficiency of foreground 
subtraction, calibration, and so on (which are expected to be 
more difficult at lower frequencies). We have not accounted 
for these effects in our forecasts, so the hgures for MID and 
LOW should be seen as a ‘best-case’ scenario. 

Note that some of the conhgurations have the same sensi¬ 
tivity at a given redshift, but exhibit small differences in their 
H (z) constraints (e.g. compare the Updated and Alternative 
conhgurations of SKA 1-LOW in Fig. 4). This is caused by 
correlations between H{z) and (Tnl. the nuisance parameter 
representing non-linear effects. The SKAl-LOW Alternative 
conhguration constrains ctnl better due to its extended red¬ 
shift coverage, for example, which results in a slight improve¬ 
ment of the marginalised H{z) constraints. 

4.2. Dark energy equation of state 

Fig. 5 shows corresponding forecasts for the equation of 
state parameters wq and Wa, including Planck CMB and 
BOSS ESS Fisher matrix priors. The curvature, Vln, and 
growth index, 70 , have been marginalised over. A summary 
of the ID marginal constraints on wq and Wa is also given in 
Table 4, assuming different sets of growth parameters. 

The sensitivity of SKAl-LOW to H{z) and Da{z) is rela¬ 
tively low, so it has little to add over BOSS in terms of con¬ 
straints on the equation of state parameters. SKAl-MID Band 
1 (Alternative conhguration) offers a signihcant improvement 
thanks to its stronger (sub- 1 %) constraints on the background 
expansion and more suitable redshift range - leading to a 
5 — 6 % constraint on u>o - but as anticipated by Fig. 4, the 
Ha survey outperforms it, reaching 3 — 4% precision on wq. 
The difference is not as drastic as one might expect, however, 
mostly thanks to MlD’s lower Zj^in compensating for its lower 
sensitivity at higher z. As expected, the SKA2 HI galaxy sur¬ 
vey consistently offers the tightest constraints (a(wo) « 0 . 02 , 
a(wa) « 0.08) because of its high sensitivity, large survey 
volume, and essentially ideal redshift range. 

Most interesting is the SKAl-MID Band 2 IM survey, 
which yields exceptionally strong constraints on (wq, Wa), ap¬ 
proaching those of the SKA2 HI galaxy survey, or even sur¬ 
passing it for some parametrisations. This is mostly due to its 
ability to break certain degeneracies by reaching the very low- 
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Fig. 6 .— Left panel: Forecast constraints on f<j% for the vaiious SKA galaxy surveys (GS) and intensity mapping surveys (IM). For SKAl-MID surveys, 
dashed lines show the ‘Alternative’ configuration. Unfilled circles show the errors on existing measurements of /erg from large scale structure surveys, taken 
from the compilation in Macaulay et al. (2013). Right panel: Correlation coefficient, r(x, y) = (Jxy I (crx<^y)-, between the growth rate {fers) and bias (berg) as 
a function of redshift. 


est redshifts. While the forecast errors on H{z) and Da{z) 
are larger than those of the SKA2 galaxy survey by a factor 
of 2 — 3, the high precision of MID Band 2 at z « 0 yields 
tighter constraints on other parameters, notably Hq, that wq 
and Wa are correlated with. Adding a tight (but achievable; 
Suyu et al. 2012) external prior on Hq at the ^ 0.5% level 
breaks the degeneracy for other surveys too, reducing (but not 
completely removing) SKAl-MID’s advantage over the Ha 
galaxy survey for example. 

Note that the constraints on wq and Wa are insensitive to the 
choice of growth parametrisation. Even the relatively loose 
constraint on Wa from SKAl-MID Band 1 changes by less 
than one percentage point across all four growth parametrisa- 
tions in Table 4. 

Now that we have an overview of how well various exper¬ 
iments can measure the dark energy equation of state, it is 
pertinent to ask how effective these constraints will be in actu¬ 
ally testing realistic dark energy and modified gravity models. 
This is a difficult question to answer however, as the space 
of plausible models is extremely broad and complex (e.g. 
see Clifton et al. 2012), and many models allow expansion 
histories that are arbitrarily close to GRh-ACDM. Taking the 
simplest class of minimally coupled scalar field quintessence 
models. Marsh et al. (2014) showed that there is no ‘target un¬ 
certainty’ on Wq and Wa beyond which most alternative mod¬ 
els could be distinguished from each other or ruled out. There 
is therefore no particular precision level to aim for - we can 
only continue to improve w{z) constraints in the hope that 
some deviation from GRh-ACDM will turn up to give us more 
clues about what to look for. The role of the SKA will be to 
substantially improve constraints on w{z) in the key low red- 
shift regime, extend measurements to higher redshifts with 
LOW, and then to provide the highest precision observations 
at z < 2 with an SKA2 HI galaxy survey. 

4.3. Linear growth rate 

Forecasts for the growth rate, fag,, are shown in Fig. 6. Ex¬ 
isting constraints from the literature are also shown, based on 
the compilation in Macaulay et al. (2013). Conservatively, the 
bias has been marginalised as a free parameter in each redshift 
bin in all of the forecasts, as described in Sect. 3.6. 


The existing constraints are restricted to the lowest redshifts 
only, z < 0.8, and are at the 10 - 20% level. The SKAl-MID 
galaxy surveys mildly improve on this, reaching the few per¬ 
cent level at the same redshifts. The MID Band 2 IM survey 
is again much more powerful, yielding 0.5 — 0.6% constraints 
at z > 0.2 that are insensitive to the assumed configuration. 
A precision of 0.6 — 3% is achievable with the MID Band 1 
IM survey over 0.7 < z < 2, but again this is substantially 
bettered by the Ha survey, which reaches 0.4% over much of 
the same redshift range. 

The SKA2 galaxy survey again provides the tightest con¬ 
straints over a wide redshift range, reaching the ^ 0.3% level, 
although the Ha survey can slightly surpass it at z « 2. At 
z > 2, HETDEX can place ^ 4% constraints on fag out to 
z « 3.5, which is roughly a factor of 2 better than the Alter¬ 
native configuration of SKA 1-LOW. Depending on the final 
band specification, LOW can put ^ 10% constraints on fag 
over the entire redshift range 2 < z < 6 however, which is 
beyond the capabilities of any other survey. 

As mentioned previously, the bias is an important source of 
uncertainty in large-scale structure analyses. The right panel 
of Fig. 6 shows the correlation coefficient between fag and 
the bias, bag, in each redshift bin. In all cases the correla¬ 
tion is moderate, except for the SKAl-MID Band 2 IM sur¬ 
vey at low redshift, where a strong anti-correlation arises. We 
have already been conservative in marginalising over the bias 
in each redshift bin, so we do not expect our predictions to 
strongly depend on the assumed bias model. 

As with the equation of state, it is useful to consider how 
well the fag constraints will be able to distinguish between 
different dark energy and modified gravity models. Perenon 
et al. (2015) randomly generated thousands of Horndeski EFT 
(effective field theory) models subject to the condition that 
they meet a set of viability criteria, and then calculated the 
distribution of fag{z) functions that the models predict (see 
their Fig. 4). While their analysis is not entirely general (it 
depends on a particular parametrisation for the evolution of 
the EFT coupling functions), it does give some idea of the 
‘typical’ range of fag{z) for a broad class of MG models, 
and so we will use their results for illustration. 

According to Perenon et al. (2015), for z « 1—2 a precision 
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TABLE 4 

Marginal errors (la x 100) on cosmological parameters for various parametrisations, calculated for four experiments: 
SKAl-MID B1 Alt. (IM); SKAl-MID B2 ALT. (IM); SKA2 (GS); AND A H« SURVEY (GS). 


Fiducial = 

0.005 Mpc-i 

0.01 Mpc-i 

0.02 Mpc-i 

Survey 


log 


log kfj. 


log 

LOW Alt. 

0.430 
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0.430 
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0.03 

MID B2 Rebase. 

0.094 

2.27 

0.095 

1.57 
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Fig. 7.— Constraints on fag for several surveys, binned by scale. 

of ~10% on fas would be required to start ruling out a signifi¬ 
cant number of models from the general Horndeski 4/5 class^ 
(see that paper for a definition). This target is well within 
the reach of the SKAl-MID Band 1 IM survey, and Ha and 
SKA2 galaxy surveys. The more restrictive Brans-Dicke sub¬ 
class predicts much less variation in fag in this range though, 
requiring ~ 0.5% precision to make any headway. The picture 
at low redshift is more encouraging however, with a ^ 5% 
constraint at z = 0.4 being sufficient to rule out a sizeable 
fraction of models from both the Brans-Dicke and Horndeski 
4/5 classes. The models predict a large spread of values of 
fas at z < 0.2, so the SKAl-MID Band 2 IM survey could 
be particularly powerful in constraining these theories. 

Finally, we present constraints on the growth rate as a func- 

^ This class includes all Horndeski scalar field theories, including Brans- 
Dicke, f{R) gravity, quintessence, K-essence, Galileons and others. 


TABLE 5 

la ERRORS ON Afj, AND logj^Q FOR A FIDUCIAL AMPLITUDE 
= 0.01 AND SEVERAL HDUCIAL VALUES OF ONLY THE PLANCK 
PRIORS ARE INCLUDED. 

tion of scale in Fig. 7. The bias, has, was also binned in 
scale and marginalised over. On linear, quasi-static scales 
(0.01 < fc < 0.1 Mpc“^), the constraints from all four sur¬ 
veys that we consider (SKAl-MID Band 1 and 2 IM surveys, 
and the Ha and SKA2 galaxy surveys) are broadly similar - 
all at the 1 — 2% level for z « 1. This is in contrast to the 
scale-independent constraints, where the differences between 
surveys are more pronounced. This is because of the better 
sensitivity of the galaxy surveys to smaller scales (k > 0.1 
Mpc“^; not shown), where there are more Fourier modes - 
autocorrelation-mode IM surveys are less sensitive to small 
scales due to their limited angular resolution. On large scales 
(k < 0.01 Mpc“^) the constraints are worse by an order of 
magnitude (~ 30%, instead of ~ 1%), but otherwise follow 
a similar pattern. Note that MID Band 1 does become more 
sensitive than the Ha survey at z « 2 on these scales due to its 
larger survey area, although the foreground removal required 
by IM experiments is least efficient at small k, potentially re¬ 
moving this advantage (galaxy surveys also suffer from large- 
scale systematics though; see Gil-Marin et al. 2015). 
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4.4. Growth parametrisations 

Table 4 shows marginal constraints on several parameters, 
for four different growth parametrisations. These can be com¬ 
pared with previous works (e.g. Di Porto & Amendola 2008). 

As with the equation of state parameters, an SKAl-MID 
Band 2 IM survey is capable of putting remarkably strin¬ 
gent constraints on the growth parameters, rivalling the con¬ 
siderably more precise SKA2 HI galaxy survey for all four 
parametrisations. This is again an artefact of its high precision 
at the lowest redshifts, which strongly breaks degeneracies 
with parameters like Hq. Adding the BOSS prior reduces this 
advantage slightly, bringing the Ha survey results to within a 
factor of ^ 2 of the MID Band 2 constraints. 

Forecasts for the ‘theoretically-motivated’ parameters 
(A^, k^) are presented in Table 5 for a small hducial devi¬ 
ation amplitude = 10“^ (i.e. a 1% correction to the Pois¬ 
son equation at z = 0), and several different scale parameters, 
kfj,. The Planck prior has been included (but not the BOSS 
one), and wg, Wa, and have been marginalised as usual. 
The two MG parameters are correlated; hxing typically 
improves the forecast errors on by a factor of a few. 

None of the surveys will be capable of detecting a deviation 
from GR at this level on their own unless is quite large, 
suggesting that a multi-tracer approach (McDonald & Seljak 
2009) may be necessary to improve precision. Even the SKA2 
constraint of (t(A^) = 0.017 for k^ = 0.01 Mpc”^ falls 
shy of a lu detection. Fixing wg, Wa, and fix also makes 
little difference. This is because the error on A^ is mostly 
driven by the bias uncertainty - hxing the bias improves the 
SKA2 + Planck constraint by roughly a factor of two. For 
comparison, a full-sky cosmic variance-limited galaxy survey 
from z = 0—3 (with b{z) = Vl + z) would achieve (t(A^) = 
0.004 with or without the bias marginalised. We conclude 
that, while a multi-tracer analysis should be helpful, the more 
immediate limitation is probably the bias rather than cosmic 
variance. 

In Simpson et al. (2013), the RSD-only constraint on their 
p,o parameter (similar to our A^) is based on hve ~ 20% mea¬ 
surements of fas at z < 0.8 from WiggleZ and 6dFGS. They 
hnd a^o = 0.25 for a ACDM background. Planck Collab¬ 
oration (2015b) improve this to =-0 21 adding more 
datasets (see their Table 6, for the parameter /tq — 1). For com¬ 
parison, our BOSS forecast (based on the BOSS specihca- 
tion in Font-Ribera et al. 2014) predicts a{A^) = 0.04 (with 
hxed k^ = 0.01 Mpc“^) from seven redshift bins with ~ 4% 
measurements of /erg. The actual BOSS dataset yielded RSD 
measurements in only a couple of redshift bins however, with 
^ 10% errors on fas (Tojeiro et al. 2012; Samushia et al. 
2014). This gives some idea of the level of optimism of our 
forecasts; real analyses can lose a great deal of sensitivity 
compared to the ‘ideal’ Fisher approach, e.g. due to cuts and 
additional modelling uncertainties introduced as a by-product 
of correcting for systematic effects. While one can hope that 
analysis methods become less ‘lossy’ in future, it is likely that 
systematics will in fact get more complex and numerous as 
raw observational sensitivity improves. Our forecast preci¬ 
sion levels should therefore be considered achievable only up 
to a factor of a few. 

Finally, the scale dependence can be constrained to within 
an order of magnitude by both SKA2 and the Ha survey, re¬ 
gardless of the hducial scale k^. The Band 1 SKAl-MID IM 
survey is better at measuring this parameter than Band 2 be¬ 
cause of its wider (and higher) redshift range, meaning that 


it can access a larger volume and thus sample large scales 
better. The constraints degrade signihcantly for the smallest 
hducial as the scale-dependent modiheation is shifted to 
larger scales where fewer modes can be measured. It was 
argued in Baker & Bull (2015) that a much larger scale of 
~ "H « 2 X 10“^^ Mpc“^ is the most natural choice in 
the absence of a specihe MG theory, so the hducial values 
we have chosen here are not particularly well-motivated. A 
proper treatment of scale-dependent MG effects on ultra-large 
scales requires the relaxation of several assumptions we have 
made in our forecasting formalism, and so we will leave this 
to future work. 

5. DISCUSSION 

Forthcoming large-scale structure surveys will greatly ex¬ 
tend the range of scales and redshifts over which General Rel¬ 
ativity can be tested. In this paper, we examined how surveys 
of neutral hydrogen with the Square Kilometre Array (SKA) 
- using either spectroscopic galaxy redshifts or 21cm inten¬ 
sity maps - can be used to probe dark energy and modihed 
gravity theories from z = 0 to z « 6. 

Using galaxy clustering and redshift-space distortions as 
our observables, we performed Fisher forecasts for the ex¬ 
pansion rate, angular diameter distance, and growth rate as a 
function of redshift and scale, which we then mapped onto 
several phenomenological parametrisations of modihed grav¬ 
ity and dark energy theories. Two possible designs - ‘Re¬ 
baselined’ and ‘Alternative’ - were considered to account for 
uncertainty in the SKA Phase 1 speciheations. 

By way of an ‘executive summary’, we will now discuss the 
implications of our forecasts for each of the proposed SKA 
cosmology surveys in turn; 

• We conhrm the hnding of Bull et al. (2015a); Raccanelli 
et al. (2015) that an SKAl HI galaxy survey will not sig¬ 
nihcantly improve existing low-z constraints from optical 
surveys on its own - its maximum redshift and mean num¬ 
ber density are simply too low. It may be useful if cross- 
correlated with other surveys though; a multi-tracer anal¬ 
ysis can lead to large gains in precision of the growth rate 
measurement if the ratio of the biases of the cross-correlated 
tracer populations is large (McDonald & Seljak 2009). The 
low bias of the HI galaxies {b « 0.7 — 0.9) may there¬ 
fore be a useful resource. Multi-tracer analyses also require 
high number densities (achievable for SKAl-MID only at 
z < 0.3) and substantial survey area overlap (possible for 
most future LSS and CMB surveys thanks to the SKA’s lo¬ 
cation in the southern hemisphere). 

• A more promising prospect at low redshift is an SKAl-MID 
Band 2 IM (intensity mapping) survey, which could map 
out the large scale HI distribution over 25,000 deg^ out to 
z « 0.6 — 0.8 with high signal-to-noise in 10,000 hours. 
This would require the MID dish array to operate in auto¬ 
correlation mode, which has yet to be tested on a large ra¬ 
dio array, and which will require a considerable amount of 
further development work to deliver sufficiently good cal¬ 
ibration and noise properties (Bull et al. 2015b). Though 
more risky from this perspective, a MID Band 2 IM sur¬ 
vey offers the best prospect of delivering cosmological con¬ 
straints with Phase 1 of the SKA that are competitive with 
contemporary experiments like Euclid; our forecasts show 
that it outperforms the Ha survey for almost all of the dark 
energy/modihed gravity parameters we considered (see Ta¬ 
ble 4). This is mainly due to its high precision at z « 0, 
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which helps to break degeneracies with parameters like Hq 
- suggesting that further improvement of low-z constraints 
is very useful for testing alternative theories in general. 

• An SKAl-MID Band 1IM survey appears to be less com¬ 
petitive, reporting constraints on the growth and expansion 
rates that are a factor of a few worse than a Ha spectro¬ 
scopic survey over the same redshift range. Still, the higher 
Zmax, large survey area, and therefore substantially larger 
survey volume of the MID Band 1 survey may make it an 
important part of a future ‘multi-tracer’ strategy to detect 
novel relativistic effects or the signatures of primordial non- 
Gaussianity in the clustering of matter on ultra-large scales 
(Alonso & Ferreira 2015; Fonseca et al. 2015), or scale- 
dependent modifications to GR (Baker & Bull 2015). 

• An interferometric IM survey with SKAl-LOW would 
do little to improve constraints on the modified gravity 
parametrisations that we have adopted, mostly because it 
only covers z > 2. It will also produce less competi¬ 
tive growth and expansion constraints than HETDEX in the 
2 ^ 2 : < 3.5 range. Its main use would therefore be to ex¬ 
tend measurements of fas{z) and H{z) to 3.5 ^ ^ ^ 6, 
which may be useful for constraining theories with early 
dark energy evolution, or coupled dark energy/dark matter. 

• While concrete specifications are yet to be provided, an 
SKA Phase 2 HI galaxy survey looks set to become the 
‘best in class’ Stage IV spectroscopic galaxy survey (Al¬ 
brecht et al. 2006), thanks to its high predicted number den¬ 
sity from z « 0 to z « 1.6. Due to its broader redshift 
range and wide (30,000 deg^) survey area, the SKA2 sur¬ 
vey outperforms the Ha survey (and all other surveys that 
we considered) in practically every measure, but as shown 
in Figs. 4 and 6 the improvement is not particularly dra¬ 
matic - gaining at most a factor of two over the Ha survey 
forecasts for H{z), Dj^{z), and fas{z) where their redshift 
coverage overlaps. This does translate to substantially bet¬ 
ter constraints on the dark energy and modified gravity pa¬ 
rameters however (Table 4), typically by a factor of 2 — 3, 
although more dramatic gains are likely to require a multi¬ 
tracer approach. 

Importantly, our forecasts establish that the cosmological 
survey performance of SKAl should not be strongly degraded 
following the ‘rebaselining’ described in McPherson (2015). 
The differences between the results for the Rebaselined and 
Alternative configurations in Figs. 4 and 6 are also small - al¬ 
though the consequences of altering the frequency bands are 


potentially more significant, as the redshift range of the sur¬ 
veys is clearly an important factor in the constraints on the 
modified gravity/dark energy parameters. 

In all of our forecasts, we assumed a total survey time 
of 10,000 hours (approximately 14 months of usable on-sky 
time), with survey areas chosen to give more or less optimal 
constraints on the dark energy equation of state (see Bull et al. 
(2015b) and Yahya et al. (2015) for optimisations of the IM 
and HI galaxy surveys respectively). This is very large allot¬ 
ment of time for a general-purpose radio observatory such as 
the SKA, as many other science goals besides cosmology will 
be competing for the available observing time. While it is 
likely that a substantial fraction of the first five years of oper¬ 
ation of SKAl will be dedicated to large survey programmes, 
it may be necessary to perform the bigger surveys ‘commen- 
sally’ (i.e. at the same time). This will inevitably introduce 
some tension into the choice of survey parameters, with some 
science goals likely preferring shallow, wide surveys, and oth¬ 
ers requesting deeper, narrower surveys. Reduced time allo¬ 
cations will also impact on performance. To give some idea of 
how different survey design assumptions would affect our re¬ 
sults, and to inform commensal survey designs, we have pre¬ 
sented some measures of survey performance as a function of 
allotted time and area for SKAl in Appendix A. 

Finally, note that in this paper we considered only a limited 
subset of possible cosmological tests of GR and dark energy: 
measurements of the background expansion, using distance 
indicators like BAOs and the broadband shape of the power 
spectrum; and measurements of the growth rate of perturba¬ 
tions on linear, quasi-static scales, using redshift-space dis¬ 
tortions. There are a number of other promising tests that the 
SKA may be able to perform however, such as: low-redshift 
Tully-Fisher peculiar velocity surveys with HI galaxies (Koda 
et al. 2014); radio weak lensing*^ surveys (Brown et al. 2015; 
Harrison cfe Brown 2015; Pourtsidou et al. 2015); observations 
of relativistic effects on ultra-large scales (Lombriser et al. 
2013; Alonso et al. 2015; Alonso & Ferreira 2015; Baker & 
Bull 2015; Fonseca et al. 2015); ‘real-time’ redshift drift mea¬ 
surements (Kloeckner et al. 2015); and strong field and grav¬ 
itational wave tests using pulsars (Shao et al. 2015). 
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APPENDIX 

A. DEPENDENCE ON SURVEY PARAMETERS 

In this appendix we briefly discuss how the performance of the SKAl surveys depends on area, S'area, and survey time, ftot- 

First we consider the SKAl HI galaxy survey. Fig. 8 shows the maximum redshift, Zmax, at which the non-linear scale is 
sample variance-limited (i.e. where the signal to noise drops to unity at so that n(zinax) b^{Zjnax)P{kNLt Zmax) = 1) as 

a function of survey time and area. Beyond this redshift, shot noise begins to dominate on progressively larger linear scales. The 
comoving volume corresponding to this redshift is also shown in the right panel of Fig. 8. Note that z^ax —> 0 in the low ftot, 
high S'area region. 

For a given observing time, one typically wants to maximise the survey volume over a given redshift range, in order to reduce 
sample variance. While larger survey areas give higher volumes for a given ftot for SKAl, the maximum redshift can be relatively 


^ Weak lensing is particularly complementaiy to RSDs, as it provides 
orthogonal constraints in the {fi, S) plane, where E is the modification to 


the lensing potential, + ’1^ (Amendola et al. 2013b; Simpson et al. 2013; 
Leonard et al. 2015). 
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Fig. 8 .— The maximum sui'vey redshift and coiTesponding comoving survey volume, Vsur(0 < 2 < ^max), in Gpc^, as a function of survey time and survey 
area, for an SKAl-MID Band 2 HI galaxy survey in the ‘Alternative’ configuration. The fiducial survey time/area assumed in this study is shown as a red cross. 
We have assumed that the minimum frequency, can be adjusted without changing other specifications, and use a detection threshold of bcr. 




'S'area [deg^ ] [deg^ ] 

Fig. 9.— Marginal Icr error on wq for several SKAl-MID IM surveys, as a function of survey area and survey time. The Planck prior is included (but not 
BOSS), and 70 is the only free modified gravity parameter. 


small, which is problematic if the aim is to measure H, Da, and fag as a function of redshift. Our fiducial choice of 5,000 deg^ 
attempts to strike a balance, obtaining a relatively large volume while maintaining a reasonable Zmax of 0.3. 

Fig. 9 shows the forecast marginal Icr error on wq for the SKAl-MID intensity mapping surveys, as a function of area for 
several survey times. While only wq is plotted here, qualitatively similar results are obtained for other parameters (e.g. 70 ). 
Apart from for the shortest survey time, ftot = 1,000 hrs, there is no clear optimum survey area for any of the configurations - 
the constraint on wq improves very slowly with survey area for all S'area ^ 5, 000 deg^. The dependence on survey time is also 
relatively moderate once a certain threshold is reached; for large survey areas, <j{wq) degrades by a factor of less than 2 for both 
MID B1 Alt. and MID B2 Alt. when going from ftot = 10, 000 hrs to 3,000 hrs. 


B. HI GALAXY SURVEY NUMBER DENSITIES 


In this appendix we present fits to the number density and bias functions derived for the SKA HI galaxy surveys, which are 
based on the calculations in Yahya et al. (2015). We use the following fitting functions; 


dN 

dz 


exp(-C32:) [deg ^]; 


b{z) = C 4 exp(c 5 z). 


(Bl) 


The fits were performed using a least-squares procedure, weighted by 7 dN/dz. This produced good fits at redshifts with sizeable 
number densities, but the fitting functions are poor approximations when dN/dz is small. The results are shown in Table 6 . 


Abdalla, F. B., & Rawlings, S. 2005, MNRAS, 360, 27, 
[arXiv:astro-ph/0411342] 
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Survey 

Thres. 

Cl 

C2 

C3 

C4 

C5 

-2^max 

Afgal/10® 

SKAl B2 Rebase. -I- MeerKAT 

5it 

5.450 

1.310 

14.394 

0.616 

1.017 

0.391 

3.49 


So- 

4.939 

1.027 

14.125 

0.913 

-0.153 

0.329 

2.04 

SKAl B2 Alt. + MeerKAT 

5cr 

5.431 

1.297 

14.260 

0.616 

1.018 

0.391 

3.54 


So- 

4.935 

1.024 

14.089 

0.907 

-0.089 

0.329 

2.05 

SKA2 

lOo- 

6.319 

1.736 

5.424 

0.554 

0.783 

1.084 

950 


TABLE 6 

Fitting coefficients for dN/dz and b{z) for the SKA HI galaxy surveys, for various detection thresholds (assuming Sarea = 5,000 
DEG^ AND 30,000 DEG^ FOR SKAl AND 2 RESPECTIVELY). THE MAXIMUM REDSHIFT AT WHICH THE NON-LINEAR SCALE IS SAMPLE 
VARIANCE-LIMITED (SEE APPENDIX A) AND TOTAL NUMBER OF GALAXIES INTEGRATED OVER THE FULL BAND ARE ALSO SHOWN. 
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